Umklapp scattering and electron pairing cutoff in BCS superconductors 
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In a superconductor electrons form pairs for which the end states of normal and umklapp scattering 
may overlap. This cuts electron pairing off at a phonon frequency, ujc, low compared with the 
Debye frequency, u>o- Using this insight, correct values of 2A/kBTc (average error 8.9%) for 12 
superconductive metals, including Hg and Pb, are achieved from simple BCS formalism with an 
average lo^I'^d of 0.148: Landau's idea of a Fermi liquid may cover strong-coupling superconductors. 
The cancellation between normal and umklapp scattering may be more important than a strong 
electron- phonon interaction in reaching a high critical temperature Tc. 
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Despite the enormous effort in nearly 50 years since 
BCS 01, a fundamental issue at the very heart of the 
study of superconductivity appears to have escaped at- 
tention: in a superconductor the end states of normal and 
umklapp scattering may overlap (FIG.^. The resulting 
competition for end states frustrates the electron pair- 
ing process and leads to a phonon cutoff frequency, 
which was postulated by Bardeen, Cooper and Schrieffer 
(BCS) but not specified and is often taken to be the 
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FIG. 1: A spherical Fermi surface (open circle) and a pair 
of electrons in states k and — k. This pair is scattered into 
states k -(- q + G and — k — q — G in the umklapp process, q 
being the phonon wavevector, G the reciprocal lattice vector. 
The shaded circles represent two identical spherical phonon 
zones, angles a and /3 measure the size of the sections of the 
Fermi surface being intersected by these phonon zones. We 
have a > P because a larger section of the Fermi surface is 
intersected by the phonon sphere centered at this surface (the 
other sphere is centered above the Fermi surface at height go)- 
As a result, the electron state — k — q — G must lie inside the 
phonon sphere centered at k, where all the states are involved 
in normal scattering. Furthermore, when k -I- q -I- G runs over 
the shaded phonon zone on the left, — k — q — G also runs 
over a sphere indicated in part by the broken white circle 
(centered above the Fermi surface at height qo to mirror the 
phonon sphere on the left). 
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Debye frequency, in the current literature j2|. How- 
ever Uc must be low to reflect the low onset frequency 
of umklapp scattering. From numerical solutions of the 
BCS self-consistent equation, an average uJc/ijJo of 0.148 
is required to reproduce experimental values of the en- 
ergy gap. A, for 12 superconductive metals, Hg and Pb 
included. Without considering k-dependence, the low Uc 
leads to substantially correct values of 2A/kBTc (aver- 
age error 8.9%, Tabled FIGs. HandOJ. It therefore 
appears that simple BCS theory, based on Landau's idea 
of a Fermi liquid, is applicable to at least some strong- 
coupling superconductors. Hence pairing frustration aris- 
ing from the competition between normal and umklapp 
scattering may be more important than a strong electron- 
phonon interaction in reaching a high T^. 

First we see how in a metal the onset frequency of umk- 
lapp scattering is usually much lower than lod. Consider 
the spherical Fermi surface in Fig. ^ that does not inter- 
sect the boundary of the Brillouin zone: there is a min- 
imum distance go between adjacent Fermi surfaces. We 
can prove by geometry that qo/qo — 0.02 for the sc lat- 
tice with one electron orbital per atom, or qo/qo ~ 0.28 
for the bcc lattice with two electron orbitals per atom, qo 
being the Debye momentum. If the phonon momentum 
q exceeds go then an electron in state k may be scat- 
tered into k + q to lie on the adjacent Fermi surface. It 
makes no physical difference if the reciprocal lattice vec- 
tor G (parallel to k for simplicity) is added to an electron 
state: we can replace k -I- q with k -|- q -I- G, which lies 
on the original Fermi surface, and is known as the end 
state of umklapp scattering; it is usually far away from 
the initial state of the electron. 

Umklapp scattering often plays a significant role in 
solids at low temperatures Q , and superconductors may 
be no exception. Let k and — k be the states of a 
Cooper pair, which by nature are on opposite sides of 
the Fermi sea, like a pair of rotating doors. Correspond- 
ingly k -I- q -I- G must be paired with — k — q — G: scat- 
tering of the electron pair is synchronized, reminiscent 
of the synchronized motion of rotating doors. It is evi- 
dent from Fig. n that — k — q — G lies inside the phonon 
sphere centred at state k, where all the electron states 
are involved in normal scattering: thus the end state of 
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umklapp scattering may not be empty, a situation that 
may have serious consequences. 

In order to incorporate umklapp scattering into the 
BCS formahsm, we must make sure that in this formahsm 
normal and umklapp scattering do not drive electrons 
into the same end state at the same time. We let 
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be our ground state wave function, where is the pair 
generation operator, /ik occupation probability, |0) the 
vacuum, M and U specify paired end states of normal 
and umklapp scattering (or normal and umklapp pairs for 
short) respectively. It is easy to check that wave functions 
of the normal and umklapp pairs with the same k are 
always orthogonal to each other: the probability is zero 
for normal and umklapp scattering to share the same end 
state at the same time. The ground state energy W = 
{^\H'QG'S,\'^) is minimized if dW/dh\^ = 0, TJbcs being 
the BCS reduced Hamiltonian 1], which leads through 
the constraint = (1/2) [1 ^ ek/-B(k)] to 
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where £'(k) = [A^(k') + e|,]^, e being the energy relative 
to the Fermi level, k can be specified by either TV or U. 

According to the above self-consistent equation nor- 
mal pairs sustain superconductivity but umklapp pairs 
frustrate it. This is true only when the normal pairs 
dominate, i.e. some terms in k'(A/') survive cancellation 
by terms in k'(W). Otherwise dW/dh\^ = will be- 
come the condition to maximize W , unless we choose 
^k — (1/2) [1 + ek/£'(k)] as the constraint and conse- 
quently, we have to exchange N and U: then umklapp 
pairs sustain superconductivity but normal pairs frus- 
trate it. It appears that either normal or umklapp pairs 
alone may lead to an energy gap (/ik or 1 — ft-k interpreted 
as the occupation probability respectively), but their co- 
existence mutually cancels their effect. We are reminded 

that (1— /ik)^^'^&k — ^k^^ '^^^ vseA by BCS to generate 
excited pairs, apparently in order to forbid excited and 
ground pairs sharing the same state at the same time, 
so that the entropy of the particle ensemble could be 
counted properly 1]. No confusion need arise. If the 
normal pairs dominate, then all the terms in k'(t/) will 
be cancelled: excited pairs will be introduced in the ab- 
sence of umklapp pairs. If the umklapp pairs dominate, 

then we use (1 — /ik)^^^ — ^k^^^k generate excited pairs. 
To make the discussion concrete we assume dominance 
of normal pairs. 

A note about the Coulomb repulsion, which was 
dropped from the reduced Hamiltonian by BCS Q. Ac- 
cording to calculations based on the Bogliubov model po- 
tential, the Coulomb repulsion has negligible effect on the 



value of the energy gap in the BCS theory, compared with 
the effect of the electron-phonon interaction, if the so- 
called Coulomb cutoff frequency, introduced artificially, 
is large enough. In fact, an energy gap may arise even 
from an entirely repulsive interaction, when this inter- 
action is perturbed by an attractive interaction over a 
narrow range of phonon frequencies It appears that 
in a metal Coulomb repulsions between electrons are in a 
balance which can be toppled even by a weak attractive 
interaction. 

We may also justify the reduced Hamiltonian from the 
field-theoretic point of view. This Hamiltonian arises 
from the canonical transformation Hbcs = e~^He^, 
where H = Hc + Hp + H^.p with [H^ + Hp,S]+ H^.p = 0, 
He, Hp and Hc-p are Hamiltonians of electron, phonon 
and electron-phonon interaction, respectively, the brack- 
ets represent the operation [A, B] = AB — BA. With this 
canonical transformation we cancel H^-p in the first order 
expansion of i?BCSi4l- Apparently we can do the same 
when H = Hc+Hp+Hcoi+He-p (i.e. to use He+Hp+Hcoi 
to replace H^+Hp in the above canonical transformation, 
Hcoi being the Hamiltonian of the Coulomb repulsion), 
provided that we use the eigen- functions of Hc+Hp-\-Hcoi 
(instead of the eigen- functions of He + Hp) as the base 
functions in the second quantization. Now Hbcs mea- 
sures the second order perturbation to the energy of 
the particle ensemble with the Coulomb repulsion. This 
brings little change to the BCS theory, except that in 
Vkk' the electron energy may be shifted slightly. 

Let us proceed. In principle we can identify the terms 
in k'(7V) that have survived cancellation by terms in 
k'(W), or vice versa, and go on to solve the self-consistent 
equation, a daunting task due to the incompatible sym- 
metries of the Fermi sea and the reciprocal lattice. It is 
clear from Fig. 1 that we have to find k'(A/') and k'(Z//) 
all over again whenever the direction of k changes. In or- 
der not to do so we follow BCS and assume that phonons 
are cut off abruptly at lOq- We also assume an isotropic 
energy gap function written as A(e). We will vary lOc to 
let the calculated A match its experimental value. Our 
theory will be supported if uJc/ujd << 1 holds to reflect 
the low onset phonon frequency for umklapp scattering. 

We adopt the free electron model, on the basis that the 
BCS theory is based on the principle of variation, which 
is not sensitive to errors in the trial fimction. We also 
adopt the Debye phonon model. Therefore in 



14k' 



2hLUq\Mc 
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where cjq is the phonon frequency, q = k — k' , we have 
fk' = fk + 4:ep£^x for states close to the Fermi surface, ep 
being the Fermi energy, ^ = q/ikp, x = cos6 + kp 
the Fermi wavenumber, 9 angle between k and q. Con- 
sequently we have (Swq)^ - (ck' - fk)^ = {'^^riYi^^ ^ 
x^), 5 = Swq/4eFC = (l/4C)(c^q/t^D)(SwD/eF) = 
{kF/2q){q/qD){TD/Tp) = {Z /16y/^ {To /Tp) , Td and 
Tp are Debye and Fermi temperatures. We have 5 ^ 
10~^ in Table m but in more realistic models S may not 
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TABLE I: 





2Ao" 


2A(0)'" 


5' 




2A(0)/fcflre'= 


Zn 


2.38 


2.40 


1.50 


0.167 


3.63 (3.20) 


Cd 


0.51 


1.50 


1.20 


0.245 


3.70 (3.20) 


Hg 


1.62 


16.5 


0.44 


0.124 


4.76 (4.60) 


Al 


2.74 


3.40 


1.81 


0.183 


3.54 (3.30) 


Ga 


7.92 


3.30 


1.51 


0.096 


3.51 (3.50) 


In 


0.49 


10.5 


0.62 


0.203 


3.92 (3.60) 


Tl 


0.46 


7.35 


0.48 


0.164 


3.92 (3.57) 


Sn 


2.52 


11.5 


1.25 


0.158 


3.70 (3.50) 


Pb 


1.39 


27.3 


0.60 


0.160 


4.56 (4.38) 


V 


27.9 


16.0 


1.20 


0.068 


3.95 (3.40) 


Nb 


9.80 


30.5 


1.04 


0.107 


4.25 (3.80) 


Ta 


6.72 


14.0 


0.90 


0.095 


3.95 (3.60) 


"in 


10-*eV 











''in 10-3 

'^experimental data bracketed 



be a constant. The matrix element 
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depends on the atomic potential 6V, where ip^s are elec- 
tron wave functions (spin absorbed into k and k'), M 
atomic mass, N number of atoms in unit volume, r coor- 
dinates in real space and Q the Wigner-Seitz cell 0, ■ 
Carbotte and Dynes utilized tabulated data of the Heine- 
Abarenkov pseudopotential to estimate SV and calcu- 
lated the electric resistivity p (another effect of the 
electron-phonon interaction) and energygap A (using 
the Eliashberg formalism) separately [a, Q ■ We use the 
Mott-Jones formula [5] to link SV with p and then use this 
p to express 6V: the Heine- Abarenkov pseudopotential is 
no longer needed. Letting k' in the self-consistent equa- 
tion [without terms in k'(W)] run over a phonon sphere 
centered at the Fermi surface, we find: 
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dx 



Here F is the overlap integral function j^j and 
Aq = he^npv^ /2kBTp 



(5) 



(6) 



measures the strength of electron-phonon interaction be- 
cause it is close to the value of J^w ^kk' (to a factor ~1), 
k' — k runs over the first phonon Brillouin zone, e and n 
are the electron charge and density, Tp the temperature 
(~room temperature) when p is measured, and v De- 
bye sound velocity. Apparently strong electron-phonon 
interaction arises from numerous free electrons (large n) 
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FIG. 2: A(e) curve for Sn at T = 0. The temperature curve 
of A at e = is shown in the insertion. 



scattered frequently by atoms (large p) that move quickly 
to facilitate pairing (large v). 

Wc find A(e) through iteration. The above integration 
with respect to x is of the Cauchy type . In first iter- 
ation we let A = Aq, so that the Cauchy principal value 
has an analytical expression: 
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where C is the integration constant, 
A{x) = A2 + (e + A€FS.5){e -f AeF£.x) 

n2i1/2 



[Al + {e + AeFm' [^l + i^ + ^^F^^n 



211/2 



B{x) =Al + {e- 4eFe^)(e + ^eFS,x) 

a211/2 



211/2 



+ [Al + {e-AeFm' [^l + i^ + ^^F^^?] 

which is integrated numerically. We vary S, the upper 
limit of the integration, until A(0) matches its observed 
value. The result in FIG. [21 (dashed curve) is virtually 
identical to the full numerical solution (where the Cauchy 
principal value is found numerically) in first iteration. It 
can be seen that prominent features of A(e) have already 
emerged from first iteration: further iterations improve 
the accuracy of solution but retain the physics. 

One such prominent feature is that A(e) is structured: 
it has a peak flanked by two negative dips. For Sn the 
distance between the peak and one dip is ~ 4 x lO"'^ eV 
(FIG. 121), compared with Hujd = 1-7 x 10^^ eV: the pair- 
ing effect extends fairly deep into the Fermi sea, in spite 
of the low ujc- As a result, the curve of de/dE, which rep- 
resents the tunnelling density of states, also has a struc- 
ture featuring a dip, similar to the characteristic dip in 
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FIG. 3: Theoretical and experimental values of 2A(0) /fcsTc. 

the tunnelling experiment data Indeed uj{e) is also 
structured in the numerical solutions for the Eliashberg 
equation |^[Tol|. 

Another prominent feature is that uJc/lod << 1 (0.148 
on average for the 12 metals in Table Pi. This IS some- 
what puzzling because, assuming that A/E varies slowly 
over the range of phonon frequencies, we can integrate 
the self-consistent equation over the first phonon Bril- 
louin zone and find reasonable values of A. Indeed, if we 
let A(k') w A(k) and hence £'(k') w E{k), then we have 
2E{k) « ^j^, Vkk', or A(0) Ao/2 (more or less true in 
Table ^ . However A/E is by no means slowly varying 
but has the shape of a sharp peak, which samples the 
extremely large and positive values of Vkk' , that would 
have been cancelled by the equally large but negative val- 
ues of Vkk', had A/E been slowly varying. As a result 
A(0) would exceed its experimental value vastly unless 

UJc/i^D << 1- 



In order to justify the values of lOc/lod in Table |l] 
we solve the self-consistent equation at T > 0, which 
is identical to its counterpart at T = 0, save an addi- 
tional factor tanh(£^/2fcBT) in the integrand fl|. Now 
A = A(e, T) which is a function of T also found through 
iteration [A — A(e) implies T — We let T increase in 
small steps and use converged A(e, T) to start iteration 
at the next T . We find Tc through a quadratic curve 
fit once A(0,T) < O.OIA(O). Tabled and FIG. El show 
that 2A(0)/fcBTc matches experimental data reasonably 
well (average error 8.9%). In particular the large values 
4.6 and 4.38 for Hg and Pb have been calculated fairly 
successfully. If we were to cut phonons off at lod, then 
we would have high Tc or low (on average ~2.5% of mea- 
sured values) with an erroneous 2A{Q) /UbTc (e.g. 3.63 
and 3.67 for Hg and Pb). 

In conclusion the electron pairing cutoff frequency Uc 
in BCS theory may arise naturally. Since lJc/^d << 1 
the 12 metals in Table U Hg and Pb included, can all 
be treated as weak-coupling superconductors, where the 
BCS theory can be applied to yield e.g. correct values 
of 2A{Q) /kBTc- Landau's idea of a Fermi liquid covers 
strong-couphng BCS superconductors. Strong electron- 
phonon interaction may not be necessary for high T^- Hg 
has a small Aq but large energy gap, whereas V has a 
large Aq but roughly the same gap (see Table P) . Our 
conclusions are valid for BCS superconductors. They 
may help us to understand MkB2 or other superconduc- 
tors with umklapp scattering [12, [iS Hjl ■ 
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